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Abstract 

Using quasi-time dependent semi-classical transport theory in RTA, we obtained coupled current equations 
in the presence of time varying field and based on general scattering mechanism r oc £P. We find that close to the 
Dirac point, the characteristic exponent /3 = +2 corresponds to acoustic phonon scattering, = + 1 long-range 
Coulomb scattering mechanism, ji = -1 is short-range delta potential scattering in which the conductivity is 
constant of temperature. The ft = case is ballistic limit. In the low energy dynamics of Dirac electrons in 
graphene, the effect of the time-dependent electric field is to alter just the electron charge by e — » e(l + (Sir) 2 ) 
making electronic conductivity non-linear. The effect of magnetic filed is also considered. 

I. Introduction 

Quite recently, semiconductor nanostructures have become the model of choice for investigation of electrical 
conduction. The unique two-dimensional material graphene which was first thought to be an academic material is 
not an exception. This 2D nanomaterial is fast becoming better candidate for electronic devices. Not only because 
of its noble electronic transport properties [1|, but it also promises a good future in graphene based electronics 
industry. Graphene has received a wide academic attention and serve as a bridge between condensed matter 
physics and high energy physics Q. The advantage of this planner material is that, one can easily change its 
electronic properties by introducing tunable gap in the sample |3| or changing the number of graphene planes 
|5l . It is also possible to fabricate free standing graphene sheets (4j. Intrisic superconducting states can also be 
realized in graphene [6 1. These among other things, means that the electronic properties of graphene can easily be 
tailored to fit device conditions. 

The crystal structure of graphene is made up of mono layer of carbon atoms arranged in hexagonal lattice. The 
low energy dynamics of fermions in graphene is characterized by linear dispersion, &(k) ~ vf\p\. In a cleaned 
sample, the conduction and valence bands touch at two inequivalent Dirac points located at the corners of Brillouin 
zone. To understand the low energy transport in graphene, a relevant scattering potential (mechanism) is essential. 
However, as we shall see in this report, one does not need to consider any explicit form of scattering potential. 
It was shown in [7| that Boltzmann theory with long-range coulomb scattering can account for all experimental 
findings. Especially, when the electronic density around Dirac points are normalized. Also, within the Boltzmann 
theory using random phase approximation, coulomb scattering has been predicted to be the dominant scattering 
mechanism [8 1. Several theories including Boltzmann Transport Equation (BTE) suggest a non-universal behavior 
of minimal conductivity which nonetheless coincides with experimentally observed value times, i.e n, (Ttheory = 
o~expln ll9l [T0ll . The same BTE predict other transport coefficients which agrees well with experiment lfT0l[TTl[T3ll . 

In this brief report, we reproduce transport properties of graphene. Within the BTE and with energy depen- 
dent relaxation time depending on power law, we showed dependence of these transport coefficients on frequency 
of the applied field. The remaining of this paper is organized as follows; Section II formulates BTE and pro- 
vides arguments leading to a quasi time-dependent solution (t-BTE). In section III, we used the t-BTE to derive 
coupled current equations from which we derived conductivity and other transport quantities. The conductivity 
tensor is rederived in the presence of magnetic field in IV. The last section V contains discussion, conclusion and 
recommendations . 
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II. BTE and quasi time-dependent solution 

A time time-dependent linearized BTE has the following form 11251 

^-V r T + e(E(t)-V r ( M /e)) 



d t f(p) + v(p) ■ 
where f(p) depends on t, r and p, i.e 

f(j>) = f(p,r,f). 

The group velocity, v is constant of time, Cf is the scattering term and T is the lattice temperature. The applied 
electric field has the form E(t) = EcosQt. H is the magnetic field. Exact analytical solutions of eq.Q are very 
difficult to obtain. Especially, due to the non-linearity of Cf and the fact that v can depend on time in general. In 
view of this, we adopt some approximations including relaxation time approximation where 

Cf(p) = -T(S)(f(p)-fom (2) 

r is inverse of relaxation time r. / and / are the time-dependent and time-independent Fermi-Dirac distribution 
functions. Motivated by [16| in the absence of magnetic field (B = 0), we consider a picture where the only 
time dependent quantity in eq.Q is the electric field. Note that Mensah solution considered the space term as 
perturbation. Under the above simplified assumptions together with the steady state solution [25|,|23| the quasi 
solution is 



f(P) = T f Q e~ W fo(p)dt-e£ e- T, 'v{p)- ^V r r + ( V r Qi/«j) - E{0) dt [-^j < 



(3) 



so that it can easily reduce back to |25 1 when Q — 0. 



III. Coupled currents and transport coefficients 

The sheet current for electron and energy flux in graphene are defined as [ 1 7 1 1 1 8 1 

/c= Mv££ v(/ , )/(/>) (4 ) 
p 

and 

p 

where g s , g v are spin and valley degeneracies and A is the graphene sheet area. We convert the sums in eqs.Q and 
Q to integrals following IPT71 

substituting eq.([3} in to eqs.Q, |5]l and simplifying using an energy dependent relaxation time 

t(S) = Ad 3 . (6) 

A is constant of energy with dimensions of s/fi and B is characteristic exponent which determines the specific 
type of scattering mechanism involved. The coupled current equations are 

J e = cr(Q)(- V,.<D) + S (-V r r) (7a) 

J s = TS(-V r ^) + K(-V r T) (7b) 
The measured electrochemical potential gradient is V<D = V r (fi/e) - E. The coefficients in eqs.(|7} are; 

2 

/ +I + ^-(k B T) 2 B(B+l)^- 1 +■■■ 

_ //~>\ u mm O /o\ 
CT^Q) = — — . (8) 

ul 1 - id/Vp 

Sp = ^$-{k B T)\B + l)^ + --- . (9) 
Mq 3eT 

K B = ^^(k B T)Y +l + ---. (10) 
Mq 3eT 
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o- mi „ is the minimal conductivity in graphene expressed as cr min = 2e 2 /h and uq is a constant given by = h/A 
and has energy dimension. 

Now, to derive a particular type of scattering mechanism, we consider specific cases when r oc £ 2 , t oc £, 
t oc and t ~ constant which correspond to/? = +2, +1, -1,0 respectively. Because the electronic conductivity 
cr is the only coefficient depending on frequency, we look at its various forms for the specified f3 values. For 
/3 = — 1; the conductivity in eq.(|8]l assumes the form 

cr_i(Q) = ^- (11) 

u\ 1 + (OAyU) 2 

This conductivity is characterized by short-range potential that has the form contact (or delta) potential and may 
be due to localized impurity (defects) in the sample iflOl [HI [T9l [T2l . For /3 = 0; we get 

a (n)=^f - -g- (12) 

which corresponds to coherent ifTUl or random Dirac mass lfl2l scattering, and describes the ballistic limit for 
electronic conductivity IfTTl in graphene. Finally, for/3 = +1; the electronic conductivity becomes 

i? + ^-(k B T) 2 

cr +1 (Q)=^ 3 (13) 
m 2 1 + (QAyu) 2 

This is an important and dominant scattering mechanism in graphene [20|. It is characterized by unscreened long- 
range Coulomb (charged impurity) scattering [7|. The second term in eq.([T3|> is inevitable at finite temperatures. 
This extra term was missing in IflOl . It is the contribution due to scattering by phonons. It is clear from eq.([T3]l 
that conductivity departs slightly from linearity behavior. Fig[T]depicts this situation. 
The conductivity for acoustic phonon scattering is obtained at /3 — +2, 

frv , Cmm ^ + n 2 (k B T) 2 fi 

o"+2(^) = —j : , 7 ■ (14) 

m 2 1 + (QA/i) 2 

A. Resistivity, thermal conductivity and thermopower 

In this section we turn to eqs.|7]i to compute other transport properties of graphene. Specifically, we will calculate 
the resistivity p, thermal conductivity, k and thermopower, So for /3 = +1. We will drop the f3 subscript in the 
following equations. By inverting eqs.|7]), one can find these quantities that experimentalist usually like working 
with. We will write our equations similar to the format in [ 10]. The resistivity is, 

u\ 1+(QA^) 2 

p(Q) = j ' ( 15 ) 

O'min , 7T , 

A< + y (k B TY 



thermal conductivity 



2 n l 2 ^ 8 n 4 4 1 + (Q.A/J.) 1 



k(Q.) = -—(K B TY^ - — —(K B T) ; ,SL, (16) 

1 + 



3W X """' m 2 9hT y " < n 2 {K B Tf 



V 2 



and thermoelectric power 

c r ~ 2 n 2 (Ki.7-) 2 l + (OA^) 2 
5(,(Q) "3^ ^ l + xHK B Tf 



. ^ 2 

The new physics emerging from these equations is the linear dependence of these quantities on Q. 2 . Note that 
electron density n dependence in our equations is self manifest, since one can easily incoperate it through IfTUl , 
|[T2l n oc fi 2 at zero temperature or 

' ' tf + ^KiTf) (18) 



n (hvp) 2 
for finite temperatures. 
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Figure 1: Normalized conductivity is plotted against 0/£2 c at fixed values of doping; /i = 
0.8eV, 0.5eV, 0.2eV and O.leV respectively. 



IV. Magnetotransport 

The BTE for non-zero magnetic field is realized from eq.Q by making the transformation E — > E + v X H or 
adding the term (v x H) ■ Vg(p) in the linearized BTE. This simple replacement will not yield a general solution, 
because of the presence of the cross product. It ensures that v-(vxH) = 0. To find the general solution, one usually 
obtains separate solutions for magnetic and electric fields and superimpose them Il24ll . Il25l . Here, we obtained the 
solution as follows; If the lattice temperature is constant of space, eq.Q becomes 

fo(p) + re'v ■ VO f-^^j = f(p) + re(v x H) ■ V p f(p). (19) 

The electrochemical potential is now defined as eVO = e'V(p/e')-e'E, with e' = e/(l -iQAp). We have assumed 
time-independent magnetic field in this section. The right hand side of eq.(|T9|> can be seen as an expansion of f(p'), 
where 

p' = p+ ^f(pxH) (20) 
\P\ 

so that 

/o + TCT.Vd)(-| =/0'). (21) 



88) 

We need to invert eq.(|2"0"|> and put it in eq.(|2T|. To do this, we make p the subject as 



p = T^fifl (p '~ T,p ' xH) ' (22) 



where 77 — rev F j\p\. Eq.pT} now becomes 



f(p) = /o + i^m p ■ m - vH x va)) (-%) (23) 
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after dropping the prime. Notice the energy dependence of rj through r, i.e rj = a&~ 1 . Where a = Aev 2 F is 
identified as the mobility in units of cm 2 V~ l s~ l . Now, to compute the electric current, eq.|23) is used in eq.Q 
with p = + 1 to get 

1 + a 2 H 2 

The presence of magnetic field vector has created off diagonal elements in the electric current density tensor. To 
compute the components of the new tensor we write eq.(|2"4|) in an indicial notation as 

cr(Q) 
1 + a 2 H 2 



J i = T^znn^ ~ <* ffyva> t ). (25) 



The longitudinal and transverse components of the magnetoconductivity tensors are 

_ 7* _ o-(Q) 
°** ~~E~ l+a 2 H 2 ' 



(26) 



J x cr(Q.)aH 

*» = % = TT^- (27) 

The rest of the components are determined through the relations <j xx = <x vy and <r xy = —cr yx . In terms of magne- 
toresistivities, eqs.p6|),([27]i are expressed as 

O-xx = 2 PXX 2 . (28) 

^-y = 2 ■ (29) 

Plx+Pxy 

with the resistivity p „ = 1/cr and the Hall resistivity p vv = aH/cr, giving 

Pxy = — 

ne 

where ne = 1/Rh and R H is the Hall coefficient. 

In general, cr is complex. For this reason, we make the replacement 1/(1 + (aH) 2 ) — > 1/(1 - zar//) such that 

1 1 + iaH 



1 - /artf 1 + (aH) 2 

The full electromagnetic conductivity for the /3 — +1 takes the form 

„2 



and 



^ 2 + ^-(^ B r) 2 

o-(Q) = ^ ? (30) 

M 2 (1 - iOA/i)(l + /or//) 



ff //(^ 2 + ^fer) 2 

o-.vv (O) = —- — — . rrN (31) 



M 2 (1 - iQAfX)(l + /orH) 

We observed the effect of crossed magnetic and electric field on graphene by plotting longitudinal conductivity 
with electric and magnetic field in Figs 2[3 



V. Discussion, Conclusion 

The advantage of this approach is that, one does not need to go through rigorous process of calculating the 
specific scattering rate T(&). For instance, unlike in [10| where finite temperature conductivity was found by 
separately calculating phonon and normal relaxation times. In this brief article, we obtained similar results without 
knowing a priori the exact form of the scattering potential. The challenge, however, is that certain material 
properties (constants), like permitivity, are not integral part of our results. Nonetheless, they can always be found 
by comparing with literature. But in this report, we do not care so much about numerical values of those constants, 
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Figure 2: Normalized conductivity is plotted against £2/£2 c at fixed values of doping; /i = 
0.8eV, 0.5eV, 0.2eV and O.leV respectively. 



we only want to demonstrate the validity and the new physics inherent in the method. In a static electric field 
Q = 0, the results obtained agrees well with what is obtained in 1 10 1, IfTTI . |[T3l , |[T4l 

Using phenomenological theory based on semi-classical BTE, we reproduce transport properties of graphene 
without knowing the type of scattering process. We found that a characteristic exponent of/? = +1 corresponds 
to charged impurity scattering and is a dominant mechanism in the absence of acoustic phonons. The graphs in 
Figs |l|2|3l may not be interesting. This is because of our approximations; the low energy behavior of carriers 
around the Dirac points. It is believed that negative differential conductivity and terahertz radiations could be 
observed when a full band spectrum obtained from tight-binding approximation is used ||2T1 . 11221 . 
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Figure 3: Normalized conductivity is plotted against magnetic field at fixed values of doping; p = 
0.8<?V, 0.5eV, 0.2eV and O.leV respectively. 
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